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THEORY  OF  MICROPOLAR  FLUIDS 

A.  Cemal  Eringen 
Purdue  University 


ABSTRACT: 

Equations  of  motion,  constitutive  equations  and  boundary  conditions 
are  derived  for  a  class  of  fluids  named  micropolar  fluids.  These 
fluids  respond  to  micro -rotational  motions  and  spin  inertia  and 
therefore,  can  support  couple  stress  and  distributed  body  couples. 
Thermodynamical  restrictions  are  studied  in  detail  and  field  equations  are 
obtained  for  the  density,  velocity  vector  and  micro-rotation  vector. 

The  system  is  solved  for  a  channel  flow  exhibiting  certain  interesting 
phenomena. 
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1.  INTRODUCTION 


The  theory  of  microfluids  introduced  by  Eringen *  deale  with 
a  class  of  fluids  which  exhibit  certain  microscopic  effects  arising 
from  the  local  structure  and  micro-motions  of  the  fluid  elements. 

These  fluids  can  support  stress  moments  and  body  moments  and  are 
influenced  by  the  spin  inertia.  The  theory  of  microfluids  are, 
however,  too  complicated  even  in  the  case  of  constitutively  linear 
theory  and  the  underlying  mathematical  problem  is  not  easily  amenable 
to  the  solution  of  non -trivial  problems  in  this  field. 

A  subclasfs  of  these  fluids  is  the  micropolar  fluids  which  exhibit 
the  micro -rotational  effects  and  micro -rotational  inertia.  This 
class  of  fluids  possesses  certain  simplicity  and  elegance  in 
their  mathematical  formulation  which  should  appeal  to  mathematicians. 
The  micropolar  fluids  can  support  couple  stress  and  body  couples  only. 
Physically  they  may  represent  adequately  the  fluids  consisting  of 
dipole  elements.  Certain  anisotropic  fluids,  e.g.  liquid  crystals 
which  are  made  up  of  dumbbell  molecules,  are  of  this  type.  In  fact, 
animal  blood  happens  to  fall  into  this  category.  Other  polymeric 
fluids  and  fluids  containing  minute  amount  additives  may  be  repre¬ 
sented  by  the  mathematical  model  underlying  micropolar  fluids. 

3  4 

Recent  experiments  with  fluids  '  containing  extremely  small 
amount  of  polymeric  additives  indicate  that  the  skin  friction  near  a 


lA.  C.  Eringen,  Int.  J.  Engng.  Sci,  £,  205  (1964). 
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A.  C.  Eringen,  "Proc.  XI  Intern.  Congress  of  Appl.  Me  eh."  Springer- 
VerlAg  (1565). 
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rigid  body  in  such  fluids  are  considerably  lower  (up  to  30 -35^) 
than  the  same  fluids  without  additives.  The  classical  liavier-Stokes 
theory  is  incapable  of  predicting  these  findings  since  it  contains 
no  mechanism  to  explain  this  new  physical  phenomena.  At  the  Naval 
Hydrodynamic  Conference  at  Bergen  last  year,  September  1964,  the 
author  suggested  that'*  the  microfluid  theory  may  contain  Just  the 
right  mechanism  required.  While  it  is  too  early  to  make  the  final 
conclusion  on  this  question,  the  problem  of  channel  flow  worked 
out  in  this  paper  is  a  positive  indication  of  this  conjecture. 

In  Arts.  2  and  3  we  give  a  resume  of  the  theory  of  microfluids 
formulated  in  Ref.  1.  The  theory  of  micropolar  fluids  is  developed 
in  Art.  4.  In  Art.  5  the  thermodynamics  of  such  fluids  are  studied 
and  the  restriction  on  the  viscosity  coefficients  are  obtained. 

In  Art.  6  we  give  the  field  equations  and  boundary  conditions  and 
present  the  similarity  parameters.  Hie  last  section  of  the  paper 
(Art.  7)  is  devoted  to  the  solution  of  the  problem  of  channel  flow 
of  micropolar  fluids. 


£ 

% 


£ 
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A 

yJ.  W.  Hoyt  and  A.  G.  Pabula,  KThe  effect  of  Additives  on  Fluid  Friction," 

U. S.  Naval  Ordnance  Test  Station  Report  (1964). 

V.  M.  Vogel  and  A.  M.  Patterson,  "An  Experimental  Investigation  of  the 
Effect  of  Additives  Injected  into  the  Boundary  Layer  of  an  Underwater 
Body,"  Pacific  Naval  Lab.  of  the  Defense  Res,  Board  of  Canada,  Ifct.  64-2. 
A.  C.  Eringen,  Proc.  3th  Symposium  on  Naval  Hydrodynamics,  Bergen, 
September  10,  1964. 
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2.  LAWS  OF  MOTION 

In  our  previous  work,  Ref.  1,  we  formulated  a  theory  of  micro - 
fluids  whose  behavior  is  governed  by  a  set  of  lavs  of  motion  and  a 
constitutive  theory.  Some  of  these  lavs  are  new  to  the  mechanics 
of  continua  and  others  are  modifications  and  extensions  of  the  well- 
known  principles  of  mechanics.  These  are 
Conservation  of  mass: 

3F  +  <P  V,k  80  in  T  (2.1) 

Balance  of  momentum: 

\i,k  *  p(f/  -  V  8  0  ln  Y  (2-2) 

Balance  of  first  stress  moments: 

*«l  -  “mi  +  \im,k  +  p(ilm  -  'atJ  "  0  lD  V  (2>3) 

Conservation  of  energy: 

*  8  *ki  vi,k +  ('u  -  tu)  vu +  \i.  w +  v* +  p*1 

in  Y  (2.k) 

Principle  of  entropy: 

q 

pr  .  on  -  '£•  >  0  in  y  (2.5) 

1 

Inequality  (2.5)  is  axioaatized  to  be  valid  for  all  independent  pro¬ 
cesses  .  In  these  equations 


► 


i 


mass  density 
velocity  vector 


-  stress  tensor 


-  body  force  per  unit  mass 

-  micro-stress  average 


the  first  stress  moments 


-  the  first  body  moments  per  unit  mass 
»  inertial  spin 

-  internal  energy  density  per  unit  mass 
■  gyration  tensor 

«  heat  vector  directed  outward  of  the  body 
*  heat  source  per  unit  mass 

-  entropy  per  unit  mass 
.  temperature 


Throughout  this  paper  we  enploy  a  rectangular  coordinate  system 
x^  ,  x2  ,  Xj  and  the  Eulerian  representation.  Fig.  1.  All  vectors 
and  tensors  are  referred  to  a  set  of  spatial  rectangular  coordinates 
so  that  no  need  arises  for  differentiating  their  covariant,  contra  variant 
and  mixed  components  from  each  other.  An  index  followed  by  a  comma 
represents  partial  differentiation  with  respect  to  space  variable  x^ 
and  a  superposed  dot  indicates  material  differentiation,  e.g. 


t,<  ■  57  '  \  *  sr  +  \,i  ’i 


(2.6) 
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Here  and  throughout  this  paper  repeated  indices  denote  summation 
over  the  range  (l,  2,  3). 

For  the  spin  inertia  we  have  the  kinematical  relation  (Ref,  1, 
eq.  5.5) 


ki 


i  ,  (v.  +  V.  V  ) 
mi  mk  nk  mn 


(2.7) 


where  i  »  1^  is  called  micro -inertia  moments  and  according  to 

the  law  of  conservation  of  micro-inertia,  they  satisfy  the  partial 
differential  equations  (Ref.  1,  eq.  2.l6) 


dl, 

-tl  +  i.  v  -  i  v 
at  km,r  r  rm  rk 


i.  v 
kr  rm 


0 


In  Y  (2.8) 


Expressions  (2.1)  to  (2.5)  and  (2.8)  are  valid  at  all  parts  of  the 
body  B  having  volume  *j/  and  surface  ,  except  at  finite 

number  of  discontinuity  surfaces,  lines  and  points.  At  the  surface 
y  of  the  body  we  have  the  boundary  conditions 


\t\  -  Xl  on  y 

\tm\  ‘  on  y 


(2.9) 

(2.10) 


where  a  la  the  exterior  norml  to  end  t^  aal  »re 

respectively  the  surface  tractions  and  surface  moments  acting  on 
We  note  that  while  equations  (2.1),  (2.2)  are  well-known  from 
the  classical  continuum  mechanics,  equations  (2*3),  (2.M  and  (2.8) 
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are  new.  The  first  two  of  these  equations  (eq.  2.5  and  2.4)  reduce 
to  classical  results^. 


p«  s  *U  vl,i  +  qk,k  +  ph 


(2.11) 


when  s  ■  0  .  Equation  (2.3)  is,  however,  much 

more  general  than  (2.11)^  and  Is  the  result  of  the  nev  principle  of 
balance  of  first  stress  moments  as  against  the  limited  axiom  of  balance 
of  moment  of  momentum  of  the  classical  theory.  Equations  (2.8)  have, 
of  course,  no  counterpart  in  the  classical  continuum  theory. 

If  we  exclude  the  heat  conduction  phenomena,  in  the  present 
theory,  the  determination  of  motion  requires  the  determination  of  the 
nineteen  unknowns. 


p(&t)  ,  ita  (£,t)  ,  vv(jc,t)  ,  vw  (jc,t) 


ki 


(2.12) 


as  agsinst  the  four  ur.kn^  ja  and  p  of  the  classical  theory. 


* 

1 

* 

| 

k 


$ 

k 


^A.  C.  Erlugen,  "nonlinear  Theory  of  Continuous  Media, "  McGraw  Bill  (1962). 
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3.  CONSTITUTIVE  EQUATIONS  OF  MICROFLUIDS 


In  Ref.  1  we  also  gave  a  set  of  constitutive  equations  for  micro - 
fluids.  For  a  non-heat  conducting  medium  these  are  expressed  as 
relations  between  (t^  ,  s^  ,  and  °frJecfrlve  quantities 


*ki 


1 

5 


+  V 


(3.1) 


VM  + 


ki 


(3.2) 


aklm  5 


Vki,m 


(3.3) 


and  p  and  i.  .  Of  these  d  is  the  rate  of  deformation  tensor 
km  * 

and  t  an(i  ft  are  two  new  tensors  respectively  called  micro - 
de formation  rate  tensor  of  second  order  and  gyration  gradient.  Both 
of  these  latter  quantities  transform  like  absolute  tensors  under  any 
rigid  motion  of  the  frame  of  reference,  i.e.  they  are  objective. 

Hence  they  are  suitable  for  use  as  the  independent  constitutive 
variables . 

For  che  present  work  we  produce  here  only  the  results  of  the 

linear  constitutive  theory  of  micro -isotropic  fluids  (i.e.  i.  ■  i  6,  ). 

Km  Km 

For  the  nonlinear  theories  the  reader  is  referred  to  Ref.  1. 
i  -  [ -ir  +  X  tr  £  +  tr  (£-£)]  I  +  2*id  +  2no(fc-£)  +  2^(fcT-£)  (jA) 


► 
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&  =  t-T  +  x  tr  4  +  n0  tr  (4-4)  u  ♦  2n  4  +  (4-4-24) 


(3.5) 


=  (7,  a  +  7„  a  +  7,  a  )  8,  ,  +  (71,  a. _  +  7.  a 

1  mrr  2  rmr  3  rrm  k l  4  trr  5 


5  rlr 


+  r6  arrP  6km  +  V  *  r8  arkr  +  79  arrk)  8ia 
+  ^10  ak/m  +  71X  akmi  +  7 12  a/km  +  71}  Bmk/  +  ^l4  aimk 


+  ^15  amik 


(3.6) 


where  £  is  the  unit  tensor  and  X  ,  ,  \i  ,  \iQ  ,  ^  H  ^  , 

and  7^  to  7^  are  the  viscosity  coefficients.  Also  tr  denotes 
trace  and  a  superscript  T  indicates  transpose,  e.g.. 


10  0 


1  =  0X0,  trbu  = 


Hie  equation  of  state  for  these  fluids  can  be  shown  to  have 


the  form 


«  =  €(n,p_1) 


(3.7) 


so  that  the  thermodynamic  pressure  tt  and  the  tesperature  0  are 
defined  by 


W  z  - 


TT*  I  ri  < 


3^  'p,i 


(3.8) 


For  a  detailed  treatment  on  thermodynamics  see  Ref.  1.  For  the  thermo¬ 
dynamics  of  micropolar  fluids  see  section  5  below. 
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For  non -heat  conducting  media,  the  nineteen  unknowns  (2.12) 
must  satisfy  the  thirteen  partial  differential  equations  obtained  by 
substituting  (2*7)  and  (3.*0  to  (3*6)  into  (2.1)  to  (2.3)  and  the 
six  equations  (2.8)  so  that  che  number  of  independent  equations  are 
equal  to  that  of  unknowns.  Equations  so  obtained  are  nonlinear  in 
the  inertia  terms  and  highly  complicated  otherwise.  The  purpose  of 
the  present  paper  is  to  give  a  new  theory  applicable  to  a  large 
class  of  fluids  falling  within  the  framework  of  the  microfluid 
theory  presented  above,  however  possessing  adequate  mathematical 
simplicity  to  make  the  engineering  problems  tractable. 
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4.  MICROPOLAR  FLUIDS 


A  microfluid  will  be  called  micropolar  if  for  all  motions 
’  “\mi  1  *  ~V4k 


Micropolar  fluids  exhibit  only  micro -rotational  effectB  and  can 

support  surface  and  body  couples.  Fluid  points  contained  in  a 

Bmall  volume  element,  in  addition  to  its  usual  rigid  motion,  can 

rotate  about  the  centroid  of  the  volume  element  in  an  average  sense 

described  by  the  gyration  tensor  v  .  No  micro -stretch  of  particles 

are,  however,  allowed  (v^  is  skew-symmetric).  Ihus  micropolar 

* 

fluids  consist  of  a  kind  of  dumbbell  molecules. 

We  now  proceed  to  show  that  a  class  of  microfluids  satisfying 
(4.1)  exists.  The  theory  of  such  fluids  is  the  subject  of  the 
remainder  of  this  paper. 

Condition  (4.1)^  implies  that 


a  .  *  -a 

kim  i km 


(4.2) 


Calculating  and  frora  (5.6)  and  equating  them  and 

using  ^.l)^  and  (4.2)  we  find  that  if  (4.1)  is  to  be  valid  for  all 
motions  we  must  have 


The  present  work  complements  our  previous  work,  Ref.  2,  on  a  similar 


subject. 
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7i  -  r2  +  74  -  r5  •  0 

7?  -  78  •  0 


r10  *  712  +  7ll  '  71J 


0 


(4.3) 


so  that 


\*m  *  (7r72)(amrAi-e/rr6km)  +  <7lO-r12,(Watal) 


+  (7l4-715)  a£mk 


(4.4) 


In  view  of  skew-symmetry  conditions  (4.1)  the  independent  number 
of  and  X^m  are  respectively  3  and  9.  Thus  it  is  natural 
to  introduce  two  new  sets  of  variables  and  by 


Vr  5  £  €rki  k t 


k  l 


=  €  ,  -  V 

rkl  r 


(4.5) 


"kr  *  "*r/m  \ira  '  \im  "  '^€fmr  mkr 


(4.6) 


where  1b  the  alternating  tensor.  Here  the  axial  vector  v 

Kim  r 

will  be  called  micro -rotation  vector  and  the  couple  stress 

tensor.  The  sign  convention  for  is  identical  to  that  of  the 

stress  tensor  and  is  shown  on  Fig.  2.  Similarly  we  introduce  micro  - 

.  * 

Inertial  rotation  and  body  couple  by 


The  couple  stress,  body  couple  and  micro-inertial  rotation  intro¬ 
duced  here  are  identical  to  those  defined  in  Ref.  6,  Art.  31. 
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0  =  -€  .  .  Ci 

r  rki  ki 


a,  =  -2e  , a 


ki 


rki  r 


i  =  -€  ,  .  i,  .  ,  i,  .  =  -2e  .  .  i 

r  ”  rki  k  l  *  ki  rki  r 


Now  multiply  (2.3)  by  €r^m  and  use  (4.5)  to  (4.8)  Since 

s  .  ■  s.  this  results  in 
m  l  fra 


(4.8) 


mrk,r  +  5kir  tir  +  p(fk  *  V  *  0 


(4.9) 


Similarly  using  (4.5)  and  (4.6)  in  (2.4)  we  may  replace  the  equation 
of  energy  by 

p«  =  \t  (vi, k  '  €Ur  Vr}  +  \l  Vt, k  +  «k,k  +  Ph  (4‘10) 

An  alternative  but  useful  form  to  (4.10)  is  obtained  by  using 

v,  .  *  d,  .  +  u>,  ,  *  d.  ,  +  w 

£,k  kl  k l  ki  kim  m 


where 


U) 


_  1 


(v,.  ,  -  V.  .  ) 


k /  -  ?  'VM  -  vi,k 


and 

the  classical  spin  tensor/  a  is  the  vert i city  vector. 


(4.11) 


P€ 


4k/  d/k  *  *k/  €kir  ^ur  +  Vr^  +  ’"k/  vl, k  +  ^k,k  + 


ph 


(4.12) 


/ 
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Hie  boundary  conditions  (2.10)  are  similarly  replaced  by 


m  .  n 
rk  r 


\ 


on 


y 


where  m^  ■  €^m  is  the  surface  couple  vector  acting 

on  .  Next  we  turn  our  attention  to  the  constitutive  equations, 
Equation  (3.4)  can  be  put  into  the  form 


\l  (-ir+Xvvr,r)eki  +  ^v(vk,l  +  vi,k)  +  %(vi,k-ck/rV1 


where  we  set 


“  *  +  5 


2 (u.  -  u  )  £  K 

1  O  V 


An  alternative  form  to  (4 .14 )  is 


k  i 


(-7 r  +  X  d  )  b,  .  +  (2u  +  k 
v  rr  kl  v  v 


^  ^ki  “  *v  €kir  ^Wr  + 


If  we  multiply  (4.4)  by  cr^m  and  use  (4.5)  and  (4.6)  this 
aquation  can  be  transformed  into 


a 


v 


v 

r,r 


v.  ,  +  7 
k,£  v 


where 

°v  !  2(712  -  V  '  Pv  !  2<72  *  V 

rv  !  2(rx  -  r2  ♦  r10  -  r12  -  7lU  +  7^) 


(U.13) 


)  (4.1*0 


(^.15) 


(**.16) 


(«.17) 


('♦.IS) 


14 


Guusnarizing  the  results:  Basic  equations  of  motion  (2.1),  (2.2), 
(4.9)  energy  (4.10)  and  the  constitutive  equations  (4.14)  and  (4.17) 
constitute  a  proof  that  the  micropolar  fluids  may  exist  as  a  subclass 
of  microfiuids  vhenever  (4.3)  is  satisfied.  The  thermodynamic  restric¬ 
tions  on  the  viscosities  are  studied  in  the  following  article 
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5.  THERMODYNAMICS  OF  MICROPOLAR  FLUIDS 

In  this  paper  we  are  primarily  concerned  with  the  non-heat 
conducting  microfluids.  In  accordance  with  the  principle  of  equi- 
presence  (Ref.  6,  Art.  44)  every  constitutive  dependent  variable 
must  be  a  function  of  the  same  list  of  variables  until  contrary  is 
shown  to  be  the  case.  In  harmony  with  this  practice  then  the 
equation  of  state  of  micropolar  fluids  must  have  the  general  form 

C  .  €(n,  p-\  du,  bki,  ak/B)  (5.1) 


The  dependence  on  i  is  dropped  since  i  «  const  along  a  material 
line.  We  proceed  to  show  that  the  dependence  of  e  on  £  ,  and 
£  can  be  eliminated  on  the  ground  of  the  second  lav  of  thermodynamics 
(2.5).  Eliminating  (q  .  +  ph)/0  between  (4.10)  and  (2.5)  we  get 


^  +  \  St  {v/,k  *  ek/r  V  +  I  "ki  v|,k  +  “5^  *  0 


(5.2) 


Using  (5.1)  this  becomes 


n  .  1  /,  1  dc*  1  de  p  .  p  /dc  •  dc  *  .  de  *  . 

pr  *  pn  (1  -  V  3T*  +  b  P  b  (5T7  k/  +  3T7  bk/  3T7"  “kin* 


dp” 


kf 


U 


kirn 


+  \\l  k  *  Sir  V  *  lrak/  Vi,k  + 


qk®,k 


*  0 


This  inequality  must  be  satisfied  for  all  independent  changes  of  n  , 


1 6 


4  >  t  >  &  and  0  v  •  SiRce  it  is  linear  in  these  quantities,  it 
cannot  be  maintained  for  all  independent  variations  of  these  quanti¬ 
ties  unless 


be 

3T 


(k*) 


be 


0  ,  p _  =  0 

^kim 


qk  =  0 


e  : 


d€  | 

3^  1  -i 

p 


pr  = 


1  de 


dkk  +  6  tkl  ^ l,  k 


€kir  Vr^  +  6  \/  V/,k  "  ° 


(5-5) 

(5.^) 

(5-5) 

(5.6) 


where  through  (2.1)  we  replaced  p  by  -pd^  •  In  (5.5)  a  parenthesis 
enclosing  indices  indicates  the  symmetric  part,  e.g. 


d€ 

ST 


uu 


_  1  /dc  ,  \ 

-  <! 


Since  any  function  €  of  a  symmetric  tensor  d^j  can  always  be 
expressed  as  a  function  d^j  w  see  from  (5.5)  that  €  must  be 
independent  of  £  ,  £  and  £  .  Using  (4,l6)  and  (4.17)  the  inequality 
(5.6)  is  further  reduced  to 

pr  2  i  [X  d  d  ,  +  (2p  ♦  k  )  d  #d  +  2k  (u.+v )(u>  +v >  )  + 

6  v  kk  tl  v'  kl  £k  v  k  k  k  k 


a  v,  ,  v.  .  +  0  v,  .  v,  .  +  y  v,  .  v.  .  ]  >0 

v  k,k  l  fl  v  k,l  1,  k  v  it  k  I ,  k  — 


(5.7) 


/ 


IT 


We  have  thus  proved 

Theorem  1.  The  necessary  and  sufficient  conditions  for  the  local 
Clausius -Duhem  inequality  (2.5)  to.  be.  satisfied  for  all  independent 
processes  are :  ( i )  e  must  he  independent  o£  4  ,  and  &  ; 

(ii)  temperature  6  and  pressure  7 r  must  be  defined  by  (5*9)  and. 

(iii)  inequality  (5*7)  must  be  satisfied  for  all  possible  motions . 
We  now  investigate  the  restrictions  emanating  from  the  satisfaction 
of  (5*7)  for  all  independent  4  ,  w  +  v  and  v.  .  .  It  is  clear 

**“  ”  by  j£ 

that  for  all  values  of  4  irrespective  of  u>  +  v  and  v  we 

by  K 

must  have  the  classical  conditions 

i  (5X  +  2m  +  k  )  >  o  ,  JL  >  o 

which  are  necessary  and  sufficient  for  the  non -negative ness  of  the 
terms  containing  4  .  Similarly  we  must  also  have 

*v/e  >  o 


in  order  that  pr  be  non-negative  for  all  values  of  w  +  v  . 
Finally  the  conditions  in  a ^  and  y ^  are  obtained  by  making 

the  last  three  terms  in  (5.7)  non-negative,  i.e. 


1 


)  >  0 


This  expression  can  be  written  as  a  quadratic  fora  in  a  nine  dimen¬ 
sional  spac*.  i.e. 
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*11  'l  °  ■  *y  ■ 

where 


yl  ■ 

V 

y 

y2 

6  v2,2 

’  y3  B  V3,3 

y4  5 

V 

y 

y5 

5  V 

'  y6  5  V2,3 

— *J 

III 

V 

y 

y8 

V3,l 

’  y9  =  vl,3 

ail  * 

a22  “ 

a33 

s 

1  K  +  ^v  + 

V  ’  al2  =  ai3 

a*5  = 

"67  = 

> 

s 

Pv/e 

a44  = 

a55  = 

a66 

= 

a77  =  a88 

a99  S  yJe 

all  other  a  =  0  . 
^  J 


Ihe  characteristic  values  ^  of  ay  are  obtained  by  solving  the 
equation 


det  (au  * a  V  =  0 


Ifte  nine  roots  for  a  are 

al  =  a2  *  a3  “  ?v  -  Pv 

-  *5  -  a.  -  a?  -  a8  .  7y  +  0y 

a9  *  +  pv  +  \ 


vA 


i 
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In  order  that  the  a  y  y  >  0  to  be  satisfied  for  all  y  we 

lj  i  j  i 

must  have 

(rv  -  Pv)/e  >o  ,  (rv  +  ev)/e  >  o 
(3av  +  Pv  +  rv)/e  >  o 

Hence 

Theorem  2.  The  necessary  and  sufficient  conditions  for  the  inequality 
(5.?)  to  be  satisfied  fur  all  motion  are 


(3\  +  2^v  +  Kv)/e  >  o  ,  nv/e  >o  ,  Ky/e  >  o 
(3av  +  2r5)/e  >  o  ,  -7v/e  <  Pv/e  <  rv/e  ,  rv/e  >  o 


(5.8) 


These  are  the  conditions  on  the  viscosity  coefficients.  In  general 
ve  also  have  6  >  0  . 

Corollary.  The  necessary  and  sufficient  condition  for  the  local 
Clausius -PubeL  inequality  to  be  satisfied  for  all  independent  processes 
are  (5.8).  This  result  is  clear  as  a  combination  of  Theorems  1  and  2. 


6.  FIELD  EQUATIONS 


The  differential  equations  satisfied  by  p  ,  v,  and  V  are  given  by 

K  IS 

(2.1)  and  combinations  of  (4.l4)and  (4.17)  with  (2.2)  and  (4.9),  i.e.. 


&  +  (pvj 


3F 


k  ,k 


(6.1) 


-TT  .  +  (A  +  U  )  V-  .  -  >  (|i  +  K  )  V.  aa  +  K  €.  -  V  . 

,k  v  v'  £,k£  Kv  v  k  ,H  v  kim  m,£ 


+  P(fk  -  vk)  =  0  (6.2) 


(a  +  6  )  V,  +  7  v.  +  K  6.  -  V  *  -  2k  V, 
v  Kv'  £,k£  'v  k,££  v  k/m  m,£  v  k 


+  p(\  -  3  \)  «  0  (6.5) 


where  a  superposed  dot  indicates  the  material  differentiation,  i.e. 


*  5T  +  Tk,f  vl 


b\ 

5t“ 


+  V  V 

k,i  £ 


(6.4) 


The  partial  differential  equations  (6.1)  to  (6.3)  are  the  field 
equations  of  the  micropolar  fluids.  Under  appropriate  initial  and 
boundary  conditions  they  are  capable  of  predicting  the  behavior  of 
such  fluids  in  a  unique  fashion.  The  existence  and  uniqueness 
theorems  must  of  course  be  proven  in  order  for  the  underlying  mathe¬ 
matical  problem  to  be  "well -posed. "  Presently  we  only  suggest  some 


21 


initial  and  boundary  conditions. 
Initial  conditions  at  t  ■  0 

p(z»0)  =  P0(x) 

vk(*’0)  = 

\U<0)  ■  vk(jc) 


where 


pQ  ,  v^  and  are  to  be  prescribed  throughout 

Boundary  conditions  at  a  rigid  boundary 


z(%t)  :  jrB 

ZUgA)  *  vB 


(6.5) 


(6.6) 


where  is  a  point  on  a  rigid  boundary  having  prescribed  velocity 

D 

y_  and  prescribed  micro -rotation  vector  v  .  Conditions  (6.6)  express 
13  ‘'"Jj 

the  assumption  of  adherence  of  the  fluid  to  the  solid  boundary. 

Boundary  conditions  involving  prescribed  forces  and  moments 
In  place  of  (6.6)  we  may  prescribe  boundary  forces  and  moments  as 
expressed  by  (2.9)  and  (4.13),  i.e. 

\t\  ■ 

(6.7) 

“ki0*  s  mi 


Other  types  of  mixed  conditions  are  possible.  The  final  Judgement 
on  these  questions  requires  theoretical  work  on  the  question  of 
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existence  and  uniqueness  and  experimental  work  on  the  flow  conditions. 

Equations  (6.1)  to  (6.3)  are  expressed  in  rectangular  coordinates. 
Vector  expressions  of  these  equations  useful  for  work  in  other  systems 
of  coordinates  are 

+  Z  '  (p  x)  =  0  (6.8) 

(*v  +  2nv  +  Ky)  2  2  •  £  -  (nv  +  *v)  2  x  2  x  i  +  «y  2  *  .X  -  2  * 

+  P t  ■  p(  -  xx  (2*  X.)  +  |s(l2)]  (6.9) 

(ay  +  Py  +  7V )  2  2  •  X  -  rv  2  x  2  x  n  +  «y  2  x  x 

-  2xy  v  +  pi  r  p  J  v  (6.10) 

•  • 
where  v  does  not  possess  as  simple  an  expression  as  x  .  There 

is,  however,  no  particular  difficulty  in  calculating  it  through  its 

tensorial  form,  cf.  [2,  17,  also  Appendix], 

We  note  that  for  Ky  «  dy  «  Py  =  =  0  and  vanishing  Jt 

through  (6.3)  we  get  vk  s  0  and  (6.2)  reduce  to  the  celebrated 

Navier-Stokes  equations .  Note  also  that  for  Ky  =  0  the  velocity 

y  and  the  micro -rotation  are  uncoupled  and  the  global  motion  is 

unaffected  by  the  micro -rotations . 

The  classical  Stokes  conditions  3A  +2 p  »  0  for  the  micro- 

v  v 

polar  fluids  have  the  corresponding  form 

3A  +  2p  +  k  *  0 
v  v  v 


(6.11) 
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to  which  we  place  no  great  faith. 

For  an  incompressible  fluids  p  =  const  ,  2  •  Z  2  0  and 

tt  is  replaced  by  an  unknown  pressure  p  to  be  determined  from 
the  boundary  conditions. 

The  similarity  parameters  of  the  micropolar  fluids  are  obtained 
by  non-dimensionalizing  equations  (6.1)  to  (6.5).  Thus  let  L  and 
T  be  respectively  some  characteristic  length  and  time  and 

Z  -  S/L  ,  t  ;  t/T  ,  £  s  v/vq  ,  v  :  v/vq 

7  •  ir/TT0  ,  p  =  p/p0  ,  f  =  f/fo  ,  J  =  J/J0  (6.12) 

are 

where  tt  ,  p  ,  v  ,  v  ,  f  and  J  /some  reference  values  of  ir  , 

O  O  O  O  O  0 

p  ,  Ixl  ,  |v|  ,  |£|  and  J  respectively.  Substituting  (6.12)  into  (6.1) 

I 

to  (6.5)  and  using  (6.4)  we  get  the  non-dimensional  equations 

n*.  —  +  (p  v  )  =  0  (6.15) 

I 

n,  v.  .  .  +  n.  v,  ..  +  n,  e,  .  v  -  n.  tt  , 

1  i,ki  2  k  ,11  5  k/m  m,£  4  ,k 

+  fk  '  n5  ^  *  7k,/  7I>  *  0  (6*l4) 

hl  v.  ,  .  +  m„  v,  . .  +  m.  €.  #  v  .  -  2m,  v, 

1  i,ki  2  k 5  kirn  m,/  4  k 

+  ?(b6  Tk  -  n5  £  ‘  V  7i>  ‘  0  (6‘15) 

1 

| 
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where 


n  z  (A  +  Li  )/p  v  L 
1  v  *v  7  o  o 


(li  +  k  )/p  vL  ,  n*  ;  k  V  /p  v  c 
Vhv  v7'  o  o  7  3  v  o/Ho  o 


n4  5  ,ro/povp2  >  n5  '•  L/Tvo  ’  n6  1  foL/vo2 


m.  E  (a  +  P  )/p  J  v  ,  m_  s  y  / p  j  v  ,  m_  =  k  / p  J  v 
l  v  v  7  o  o  o  7  2  7v7Ko  o  o  7  3  v'  o  o  c 


n».  f  k  L/p  J  v  ,  rac  =  n, 

4  -  v  /Kouo  o  7  5  3 


U/j  v  v 
o  7  o  o  o 


(6.16) 


Of  these  are  the  reciprocal  Reynold  numbers,  n^  ,  n^ 

and  n^  are  well-known  from  the  Navier-Stokes  theory.  ‘The  present 
theory  introduces  six  new  numbers  namely  n^  ,  m^  ,  ,  m^  ,  m^ 

and  m^  .  For  a  given  fluid  is  proportional  to  so  that  the 

only  new  parameters  are 

2 

n  ■  k  v/pv  m_  ■  7  /p  j  v  ,  m.  ■  k  L/p  J  v 

3  v  o'  o  o  7  2  v7  o  o  o  7  4  v  7  o  o  o 


S  =  m4/m3  '•  VoL/vo  >  “6  :  VA>vovc 


(6.17) 


Ihe  four  of  these  new  similarity  parameters  represents  the  relative 
importance  of  rotational  viscosities  to  the  inertia  terms  and  the 
fifth  m^  the  relative  micro -rotation  velocity  to  the  velocity. 
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7.  FLOW  OF  MICROPOIAR  FLUIDS  IN  A  CIRCUIAR  PIPE 

Here  ve  give  the  solution  of  the  field  equations  (6.8)  to  (6.10) 
for  a  steady  motion  of  micropolar  fluids  in  a  circular  channel.  The 
appropriate  coordinate  system  for  this  problem  is  the  cylindrical 
coordinates  (r,  0,  z)  with  z  taken  along  the  axis  of  the  pipe.  For 
a  steady  flow  ve  seek  to  determine  the  velocity  and  micro -rotation 
components 


v  b  v  =  0 
r  0 


CD  cm  m  0 

r  z 


v  »  w(r) 

(7.1) 

<P0  -  9(r) 


Equation  of  continuity  (6.8)  is  satisfied  identically  for  p  «  const. 

and  (6.9)  and  (6.10)  with  £  =  l  =  0  give  p  «  p  «  0  and 

,r  ,0 

(n  +  k  )  (rw1)'  +  k  (rv)'  «  rp  (7.2) 

V  V  f  z 

7y(v'  +  r"^v) '  -  #c  ■ w*  -  2/c^v  «  0  (7.5) 

where  a  superposed  prime  indicates  differentiation  with  respect  to  r. 
We  also  used  p  to  denote  hydrostatic  pressure  in  place  of  t r. 

From  (7.1)  we  solve  for  v* .  Hence 

-1  r  -1 

w’  "  (K,  +  *J  (**  v  *  o  P  ,)  +  Cir  (7 •**) 

V  V  V  c  fZ  1 

Next  substitute  v'  into  (7*5).  This  gives 


1 

+  -v 
r 


•-  (k2  +  * 


^)v  - 


Pr 


(7.5) 
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where 


=  (2>Jv  *  h ,  lx  j1/2 ,  P  s 


M  +  K  '  y 
V  V  V 


Kv _  dp_ 

2^v  +  *v)\.  dz" 


(7.6) 


The  general  solution  of  (7. 5)  is  found  to  be 


v  ■  A  I^kr)  +  B  K^kr)  -  Pk  2r 


(7.7) 


where  I^p)  and  K^p)  are  modified  Bessel  functions  of  first  order 
and  first  and  second  kind  respectively.  Substituting  this  into 
(7*4),  and  integrating  the  result  we  obtain 


w  -  Ky(^y  +  O'*1  k"1  [“  A  I0^kr^  +  B  K0(kr)l 
+  |  (  2\i  +  <c  )-1  p>f  r2  +  C  log  r  +  C 


(7.8) 


where  I  and  K  are  modified  Bessel  functions  of  zeroth  order  and 
o  o 

first  and  second  kind  respectively  and  C  is  an  arbitrary  constant. 

Both  w  and  v  must  be  bounded  at  r  ■  0.  Since  KQ(kr),  K^kr) 
and  log  r  become  infinite  for  r  ■  C  we  must  have  B  ■  ■  0.  We 
assume  that  the  fluid  sticks  to  the  boundary  r  ■  a,  i.e., 


w(a)  -  0  ,  v(a)  -  0  (7 -9) 

Using  (7*7)  and  (7.8)  we  determine  A  and  C  leading  to  the  solution 

»/w0  - 1  -  02  ♦  iHfV  j 
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-  P 


Il(Ao) 


(7.11) 


where 


w 

o 


2 

a 


(2py  +  Kv) 


d£ 

dz 


P  -  r/a 

X  *  k»  *  4--^— — ^ 

(i  +  K 


(7.12) 


Here  w^  is  the  maximum  velocity  in  the  classical  Poiseuille  flow 
which  occurs  at  r  *  0.  The  solution  (710)  goes  into  the  classical 
Poiseuille  flow  for  k  *  0  and  (7.11)  gives  v  *  0. 

According  to  (5*8)  with  0  >  0  we  have  and  y ^  non¬ 

negative.  Thus  ^  is  a  real  number.  For  various  values  of  7\  we 
give  on  Fig.  4  plots  of  velocity  difference  from  the  classical. 

Poiseuille  flow  and  on  Fig.  5  vh/w  .  From  Fig.  4  as  well  as  Fig. 

o 

3  we  see  that  the  velocity  profile  is  no  longer  parabolic.  Moreover 
the  velocity  here  is  smaller  than  that  of  the  classical  Navier-Stokes 
fluids.  Of  course,  micro-rotation  v  is  altogether  missing  in  the 
Navier-Stokes  theory. 

The  non-vanishing  components  of  the  stress  tensor  and  those 
of  the  couple  stress  are  obtained  through  expressing  (4.14)  and 
(4.17)  in  cylindrical  coordinates.  Hence 
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rr 


■  -  P 


t 

rz 


P 
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t  ,  a  [p.„-.&L.- 

zr  2  dz  m  +  K  I-i(X) 

V  V  «*• 

V  ■  (Pv  Vo/2a)  U  +  '  -  My  ] 

-  <VV  “re 


(7-15) 


We  note  that  t  i  t  whenever  K  i  0. 

rz  zr  v 

On  Figs.  6  to  7  are  6hown  the  surface  tractions  and  couples  on 

the  fluid  surface  adjacent  to  the  wall  at  p  =  1  for  dp/dz  <  0,  <  0 

and  of  course  7^  >  0.  The  shearing  stress  t  has  the  same  expression 

as  in  the  classical  theory.  However  the  existence  of  the  distrib- 

(Fig.  3) 

uted  couples  a  on  the  fluid  surface  /will  produce  an  effect  in  a 
thin  layer  near  the  wall,  equivalent  to  reduction  of  the  surface 
shear.  Clearly  then  the  present  theory  gives  rise  to  a  boundary 
layer  phenomena  not  present  in  the  Navier -Stokes  theory.  This  new 
boundary  layer  is  controlled  with  the  parameter  7\, 

We  believe  that  the  theory  of  raicropolar  fluids  opens  up  a 
very  worthwhile  branch  of  fluid  mechanics.  It  should  find  important 
applications  dealing  with  a  variety  of  fluids.  It  should,  in 
particular  cast  new  directions  in  the  theory  of  turbulence.  Rich 
theoretical  and  experimental  studies  are  awaiting  the  future  workers. 
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Fig.  5  Micro -Rotation 
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